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Abstract — Interference between nodes directly limits the ca- 
pacity of mobile ad hoc networks. This paper focuses on spatial 
interference cancelation with perfect channel state information 
(CSI), and analyzes the corresponding network capacity. Specif- 
ically, by using multiple antennas, zero-forcing beamforming is 
applied at each receiver for canceling the strongest interferers. 
Given spatial interference cancelation, the network transmis- 
sion capacity is analyzed in this paper, which is defined as 
the maximum transmitting node density under constraints on 
outage and the signal-to-interference-noise ratio. Assuming the 
Poisson distribution for the locations of network nodes and 
spatially i.i.d. Rayleigh fading channels, mathematical tools 
from stochastic geometry are applied for deriving scaling laws 
for transmission capacity. Specifically, for small target outage 
probability, transmission capacity is proved to increase following 
a power law, where the exponent is the inverse of the size of 
antenna array or larger depending on the pass loss exponent. As 
shown by simulations, spatial interference cancelation increases 
transmission capacity by an order of magnitude or more even if 
only one extra antenna is added to each node. 

I. Introduction 

In a mobile ad hoc network (MANET), the mutual in- 
terference between nodes limits throughput for peer-to-peer 
communication over the network. In this paper, zero-forcing 
beamforming are applied at receivers for canceling interfer- 
ence from the strongest interferers. Thereby the number of 
successful communication links per unit area, called network 
transmission capacity [1], increases significantly. 

For a multi-hop ad hoc network, the notion of transport 
capacity was introduced and analyzed in [2], which started 
a series of related studies (see e.g. [3]-[5]). Prior results on 
transport capacity typically focus on scaling laws of network 
throughput with an asymptotically large number of nodes. 
Such asymptotic results may differ significantly from the 
actual throughput of finite-size networks, and thus have limited 
practical applications. 

For single-hop ad hoc networks, the transmission capacity 
metric introduced in [6] is defined as the maximum number 
of successful communication links per unit area under signal- 
to-interference-noise ratio (SINR) and outage constraints. By 
modeling network nodes as a Poisson point process, this 
capacity measure enjoys more accurate analysis and easier 
computability compared with transport capacity [1], [6]. Trans- 
mission capacity has been used to make tractable analysis 
of various issues related to MANETs such as opportunistic 



transmission [6], successive interference cancelation (SIC) [7] 
and multi-antenna transmission [8]. 

Interference directly limits the throughput of MANETs. The 
optimal approach for reducing interference in MANETs is 
called interference alignment, which achieves the number of 
degrees of freedom equal to half of the number of interference 
links [9]. Nevertheless, this approach appears daunting because 
it requires nodes to employ jointly designed precoders and 
obtain perfect channel sate information (CSI) of interference 
channels. An alternative approach is to use physical layer 
techniques of the multiuser detection family for suppressing 
interference [10]. These algorithms, however, appear to be very 
sensitive to residual interference due to imperfect interference 
cancelation, and the near-far problem. 

This paper considers zero-forcing beamforming for inter- 
ference cancelation in multi-antenna MANETs with single- 
stream data links. Therefore, the spatial degrees of free- 
doms created by multi-antennas are dedicated for interference 
cancelation. Recently, beamforming or directional antennas 
have been integrated with the medium access control (MAC) 
protocols for MANETs to improve network spatial reuse 
efficiency (see e.g. [1 1]— [14]). Most prior work focuses on 
designing MAC protocols and relies on simulations. In [8], the 
transmission capacity for multi-antenna MANETs is analyzed, 
where interference is treated as noise and suppressed by 
averaging through beamforming. There still lacks analysis of 
the effects of spatial interference cancelation on the transmis- 
sion capacity of MANETs, which is addressed in this paper. 
Another important issue, namely how CSI inaccuracy affects 
the transmission capacity of MANETs, is investigated in a 
sperate paper [15]. 

The contributions of this paper are summarized as follows. 
This paper targets a MANET with single-stream data links 
and perfect synchronization between nodes. First, assuming 
perfect CSI, zero-forcing beamforming is applied for canceling 
interference at receivers and thereby increasing network trans- 
mission capacity. Second, based on the Poisson assumption 
on transmitting-node locations and the spatially i.i.d. Rayleigh 
fading channel model, bounds on the signal-to-interference 
ratio (SIR) outage probability are derived. Third, the scaling 
laws for transmission capacity are derived for asymptotically 
small target outage probability. Specifically, the asymptotic 
transmission capacity grows following a power law, where 



the exponent is the inverse of the antenna-array size or larger 
depending on the pass-loss exponent. 

II. Network and Channel Models 

A. Network Model 

In this paper, the locations of potential transmitting nodes 
in the mobile ad hoc network, including both active and 
inactive transmitters, are modeled as a Poisson point process 
following the common approach in the literature [1], [6]-[8]. 
Specifically, the positions of the potential transmitters form a 
homogeneous Poisson point process with the density denoted 
by A G . Each transmitter transmits with fixed probability de- 
noted by P t . Let T n denote the coordinate of the ith transmitter 
on the 2-D plane. The set $ = {T n } is also a homogeneous 
Poisson point process but with a smaller density A = PtX 
[16]. Each transmitter is associated with a receiver located at 
a fixed distance denoted as d. 

Consider a typical receiver located at the origin, denoted 
as Rq, and hence |Tb| = d. This location constraint of To 
does not compromise the generality since the transmitting node 
process $ is translation invariant. Furthermore, according to 
Slivnyak's theorem [16], other transmitters, namely $/{To}, 
remain as a homogeneous Poisson point process with the same 
node density A . 

The ad hoc network is assumed to be interference limited 
and thus noise is neglected for simplicity. Consequently, the 
reliability of data packets received by the node Ro is deter- 
mined by the SIR. Moreover, we assume that each data link in 
the network has a single stream, and communications between 
nodes are perfectly synchronized. Let S denote the random 
channel power for the link from To to Rq, and the function 
I(T n ) gives the power of interference from the T n to Rq. Thus, 
assuming uniform data transmission power for all transmitters, 

the SIR at Rq is given as SIRo = — » T ; ■ To 

simplify notation, I(T n ) is denoted as /„ in the sequel. The 
correct decoding of received data packets requires the SIR to 
exceed a threshold 9. In other words, the rate of information 
sent from a transmitter to a receiver is log 2 (l + 9) assuming 
Gaussian signaling. To support this information rate with high 
probability, the outage probability that SIRo is below 9 must 
be smaller than or equal to a given threshold < e C 1, i.e. 

Pout (A) =Pr(SIR <9) <e (1) 

where P out (A) denotes the SIR outage probability. Given e, 
P out determines the transmission capacity defined as [1] 

C(e) = (1 - e)A e (2) 

where P ou t(A e ) = e. 

B. Channel Model 

The channel model is characterized by narrow-band and fiat 
fading. Each node in the network is equipped with L antennas. 
Consequently, there exists a L x L multiple-input-multiple- 
output (MIMO) channel between every pair of nodes. Each 
MIMO channel consists of path-loss and spatially i.i.d. small 



fading components, corresponding to rich scattering. Specifi- 
cally, the channel from a node T„ to the typical receiving node 
Ro is H„ = d n a ^ 2 G n . The factor d n a ^ 2 represents path-loss, 
where d n = \T n \ and a > 2 is the path-loss exponent. The 
other factor of H n , G„, models spatially i.i.d. Rayleigh fading, 
and hence G is a L x L matrix of i.i.d. CA/"(0, 1) components. 
The above channel model simplifies analysis in this paper. 
Finally, given single-stream data links, beamforming is applied 
at each transmitter and receiver. 

III. Spatial Interference Cancelation 

Assume perfect CSI, synchronization between nodes, and 
single-stream data links. Under these assumptions, spatial in- 
terference cancelation uses zero-forcing beamforming by fol- 
lowing the procedure described in Section UlI-AI The effective 
network and channel models are presented in Section IIII-BI 
Without loss of generality, the discussion focuses on the 
typical pair of nodes To and Ro (cf. Section Hi- Al l. 

A. Algorithms 

The idea of spatial interference cancelation is to apply zero- 
forcing beamforming at Ro for canceling interference from 
strong interferers. Let f„ and vo denote the transmit beam- 
former at T n and the receive beamformer at Ro, respectively. 
From the perspective of Ro, the interference channel from T ra 
(n 7^ 0) appears as an effective channel vector h„ = H n f„. To 
facilitate our discussion, the indices of the interferers of Ro are 
sorted according to their effective interference channel norms, 
namely ||hi|| > ||h2|| > ••• > ||hf,||---. By zero-forcing 
beamforming, the beamforming vector vo of Ro is constrained 
to be in the null space of the matrix [hi,h2,-- - , hf,_i]. 
Thereby, the interference from L — 1 strongest interferers 
to Ro is canceled. Note that perfect CSI estimation of 
hi, h.2, ■ ■ ■ , ht_i by Ro is required to completely cancel the 
interference from L— 1 strongest interferers. CSI estimation at 
each receiver uses pilot symbols broadcast by transmitters. The 
issue of CSI inaccuracy is addressed in [15]. Next, an arbitrary 
transmit beamformer is applied at To, represented by fo. By 
such beamforming, multiple transmit antennas contribute no 
diversity gain. 

To avoid deep fading due to the lack of diversity gain, 
opportunistic transmission is applied as in [6]. Consequently, 
transmission at each transmitter is turned on only if the channel 
gain S is above a threshold denoted by (3, where S = 
vjHofo| 2 . It follows that the transmission probability for each 
potential transmitter is P t = Pr(5 > (3) (cf. Section Ill-Al l. 

B. Effective Channel and Network Models 

With perfect interference cancelation, Ro receives interfer- 
ence only from the nodes {T„ | n > L}. Let r n and /„ denote 
respectively the distance between T n and the origin, and the 
interference power from T n to Rq. Using this notation and 
based on the channel model in Section III-BI for n > L, 
In = PD|vjH n f n | 2 = r- a \vlG n f n \ 2 . Because both f„ 
and vo are independent of G„ and G„ is an i.i.d. complex 
Gaussian matrix, the random variable p n = |vjG„f„| 2 follows 
the exponential distribution with unit variance. Therefore, for 



n > L, the effective interference power from T n to R n is 

In P^* 

The effective power of the data link from To to Ro is given 
by S = |v H f | 2 = cT Q | vjG f 1 2 with S > (3 due to 
opportunistic transmission. Because the beamformers vo and 
fo are independent of Go as discussed in the preceding section, 
the random variable W = |vjGofo| 2 has the exponential 
distribution with the following probability density function 
(PDF) 

f w (w) = exp(-w)/P t , w > (3d a (3) 
where P t — exp(— /3cP). 

IV. Outage Probability and Transmission Capacity 
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2) Conditioned on G = g, the mean and variance of the 
secondary interference power In{g) are 



E[/n(s)] = 

Var(/ n (<?)) = 
Proof: See Appendix lAl 
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Based on auxiliary results in Section IIV-AI bounds on the 
SIR outage probability are obtained in Section IIV-BI The 
scaling law of transmission capacity is derived for the regime 
of small outage probability in Section IIV-CI 

A. Auxiliary Results 

To facilitate analysis, the interference nodes of Ro after 
perfect interference cancelation, namely {T n | n > L}, are 
separated into the strongest interferer Tl and others {T n \ 
n > L + 1}, referred to respectively as the primary and the 
secondary interferers. For convenience, denote the random 
interference power from Tl as G = II- The separation of 
interferers provides a useful result that conditioning on G, the 
secondary interferers {T n | n > L + 1} form a Poisson point 
process as shown shortly. This result is obtained by using the 
Marking Theorem [16]. To apply this theorem, a marked point 
process is defined for the secondary interferers, where the 
mark of the node T n is the corresponding interference power 
/„. Specifically, conditioning on the interference power of the 
primary interferer G = g, the desired marked point process 
is 11(g) = {{T n ,I n } | T n e $/{T },0 < I n < <?}, where $ 
is the transmitter process defined in Section III-AI Note that 
conditioning on G = g, the marks for different nodes are 
independent. Given this condition, the result in the following 
lemma directly follows from the Marking Theorem. 

Lemma 1: The process IT is a Poisson point process 
with the average number of nodes given by fi(g) = 

oo g 

2tt\ J J rp(r, dl)dr. 
o o 

This result is useful for analyzing the aggregate interference 
from the secondary interferers to Rq. Conditioned on G = g, 
this interference is I n (g) = L(T„,/„)en*( s ) I n- 

Next, the distribution of the interference power from both 
the primary (G) and the secondary interferers (In(g)) are 
characterized in the following lemma. 

Lemma 2: The interference at Rq has the following prop- 
erties. 

1) The PDF of the primary interference power G is 



B. Bounds on Outage Probability 

The analysis of the exact outage probability P ut(A) = 
E [Pr(SIRo < 0)] is infeasible due to the difficulty in deriving 
the distribution function of the secondary interference power 
In (g). Therefore, we resolve to obtaining bounds on P out 
following the approach in [1], [6], [7]. The outage probability 
can be written as 

P out (A)=E [Pr(I n (G) > Wd- x dT a -G\G,W)\ WG~ l > 
0d a ] Pr<mr 1 > 6d a ) + PriWG- 1 < 6d a ). (7) 



Thus, a lower bound of P out is given as 



Pout (A) > Pi'tW^T 1 < 6d a ). 



(8) 



This lower bound considers only the primary interference, and 
hence is tight if the primary interferer is the dominant 
source of interference. Next, an upper bound of the outage 
probability can be derived by applying the following Cheby- 
shev's inequality on ^ 



Pr(InO) > a) < min 



Var (/□((/)) 



(9) 



>-E [I n («?)]} 

Based on (0, ([8} and (0, bounds on the outage probability 
are derived as shown in the following proposition. 

Proposition 1: For perfect CSI, the bounds on the outage 
probability are given as follows. 

1) The lower bound is 

1 



PoutW 



T(L) 



E[ 1 (L,c 2 XW- 5 )] (10) 



where c 2 = 7rr(l + S)9 s d 2 , and W has the probability 
density function in (f3]). 
2) The upper bound is 



PoutW = pLW + [i - pLW] PaW 

where 

c 3 XG 2 - s 



(ID 



P Q (A)=E <^ min 



(d- 



~ X W — G — c 4 AG 2 



w 

G >'« 



(12) 



and c-i = 27rr ( A "+ 1 ) anc [ c . — 7rr ( 1 + l5 ) _ 

J a — 2 * a — 1 

The proof is straightforward and is hence omitted. The diffi- 
culty in deriving a closed-form expression for P out is mainly 
lies in that the distribution function of 7n(G), called a shot 
noise process, is unknown [6], [17]. 



C. Asymptotic Transmission Capacity 

In this section, the scaling law for transmission capacity 
is derived for small target outage probability (e — > 0). This 
scaling law also accurately characterizes transmission capacity 
in the non-asymptotic outage regime (up to 0.1) as shown by 
simulations in Section IV1 

Small target outage probability results in a network of sparse 
transmitting nodes (i.e. A — > 0). For such a sparse network, the 
useful relationship between the outage probability and node 
density is derived and shown in the following lemma. 

Lemma 3: For A — * 0, the outage probability scales with 
A as follows. 

1) For L < a, 



where Ki 



, -Pout (A) , 
Ki < lim — —f < Ki (1 

A^O X L 
_ T(l-SL,(3d c ")[TTT(&+l)0 s d 2 } 1 



K2) 



(13) 



L 

2) For 



x-2 

> a, 



Ki < lim 

A^O 



p t r(i+i) 



and K2 — 



PotltW ,. Pout (A) 

' I™o^^^ K3 (14) 



where K3 



8{c 1 d 2 8 s ) a r{-l,pd a )r(L~a+l) 
(a-2)P t T{L) 

Proof: See Appendix iBl □ 

Using Lemma |3] and the definition of transmission capacity 
in (O, the main result of this section is obtained and summa- 
rized in the following theorem. 

Theorem 1: For small target outage probability e — > 0, the 
transmission capacity scales as 
1) For L < a, 



lim 



+ n 2 )] £ ex 



> 1, 



lim ■ 

f .— n 



< 1 (15) 



where K\ and K2 are specified in Lemma \3\ 
2) For L > a, 



lim ■ 



> 1, 



lim • 



< 1 



(16) 



'El 



where K3 is given in Lemma \3\ 
The above theorem shows that as the target outage probability 
decreases, transmission capacity grows following the power 
law ae* where a and f are constants. For L > a, only bounds 
of the exponent t are known. The derivation of the exact 
exponent may require analyzing the distribution function of 
the secondary interference power (cf. Section lTH-Bb rather than 
using Chebyshev's inequality in ©. Unfortunately, no closed- 
form expression for this distribution function is known for the 
present case [17]. 

For L < a, the exponent of the transmission capacity 
power law ae* is shown in Theorem [T] to be t = 1/L, and 
a is bounded as + ^2)]^"^ < a < n 1 L . This power 

law indicates that the size of antenna array L determines the 
sensitivity of transmission capacity to the change on the outage 
constraint. To facilitate our discussion, rewrite the scaling 
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Contending node density 

Fig. 1 . Outage probability for different transmitting node densities and perfect 
CSI. The size of the antenna array is L = {2, 4}. 

law in Theorem Q] as C(e) = aei where "=" represents 
asymptotic equivalence for e — > 0. Moreover, consider two 
sets of values (C\, ei) and (C2, £2), and define the logarithmic 
ratios AC = log ^ and Ae = log Using this notation, the 
above scaling law can be simplified as 

Ae L 

The above quantity represents the sensitivity of transmis- 
sion capacity towards the change of the outage constraint. 
Its value decreases inversely with the size of antenna array. 
Specifically, computed using ( fl7l i. a hundred-time decrease 
on e reduces network transmission capacity by {10,3.2, 1.8} 
times for L = {2, 4, 8}, respectively. 

V. Simulation and Discussion 

The results in this section are obtained by simulating a 
MANET following the procedure in [18]. 

For perfect CSI, the bounds on outage probability from 
Proposition Q] and simulated values are compared in Fig. Q] 
The path loss exponent is a = 4 and the number of antenna 
per node is L = {2,4}. Two observations can be made from 
Fig. [TJ First, the bounds for L = 2 are tighter than those for 
L = 4. Second, the bounds and the simulated values of the 
outage probability converge as the transmitting node density 
A decreases. These two observations can be explained by the 
dominance of the primary interference over the secondary 
one as L increases or A decreases, where the secondary 
interference causes the looseness of the bounds on outage 
probability. 

In Fig. |2 the transmission capacity is plotted for an 
increasing number of antennas per node assuming perfect 
CSI. Furthermore, different outage constraints, namely e = 
{10 _1 , 10~ 2 , 10~ 3 }, are considered. From Fig. |2 the fol- 
lowing observations are made. First, the use of multiple 
antennas for interference cancelation leads to the increase in 
transmission capacity by an order of magnitude or more with 
respect to the case of single-antenna per node. This capacity 
gain is especially large for a small number of antennas and 
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Fig. 2. Transmission capacity by simulation for different node densities and 
perfect CSI. The size of the antenna array is L = 4 and the outage constraint 

is e = {icr\i(r 2 ,i(r 3 }. 
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Target outage probability 

Fig. 3. Comparison between asymptotic bounds on transmission capacity and 
the exact values obtained by simulation perfect CSI and the size of antenna 
array L = {2, 3, 4} 



small target outage probability. For example, for e = 10 _1 , the 
use of three antennas per node provides transmission capacity 
seven times of that for the single-antenna case. The capacity 
gain by using additional antennas diminishes rapidly as the 
number of antennas per node increases. Second, the outage 
constraint affects transmission capacity significantly for a 
small number of antennas per node. Nevertheless, transmission 
capacity becomes insensitive to the change on the outage 
constraint as the number of antennas increases. 

In Fig. |3] asymptotic bounds on transmission capacity in 
Theorem Q] are compared with the exact values obtained by 
simulation for perfect CSI and the range of target outage 
probability e 6 [1CP 5 ,0.1]. As observed from Fig. [5] the 
asymptotic upper bound on transmission capacity is very 
tight for L = {2,3} even in the non-asymptotic range e.g. 
e G [0.01,0.1]. The tightness of this bound is due to the 
dominance of primary interference for interference cancelation 
with small sizes of antenna array. For L = 4, both the 
asymptotic lower and upper bounds are tight. The tightness 
of asymptotic bounds implies the slopes of the transmission 
capacity vs. target outage probability curves are approximately 
equal to 

Appendix 

A. Proof of Lemma [2] 

Define a marked Poisson point process as Ai(g) = 
{{T n ,I n } | T n e $/{T }, I n > g}. Given that 
the node density of M (g) follows from the Marking Theorem 



p{M{g)) 



= X I / 2nrf p (t)drdt 
Jo Jo 

= TrT(S + l)\g- s , 5:=-. 

a 



The probability density function of G is obtained by differen- 
tiating the above CDF. 

Using Campbell's theorem} 16], the expressions for 
Var(In(<?)) and E[/n(.9)] are obtained as 

rOO rOO 

E[/ n | g] = 2tt\ / / r 1_Q pe~ p drdp 
Jo J(R) 1/a 



2nXT{S + l) 1-4 
a-2 9 ' 



(18) 



Var(Jnlff) = 2ttA / / r (r~ Q pf ' 'dpdr 

Jo J(R) 1/a 



(f) 

7rr(3 + l)A g2 _^ 
a — 1 



(19) 



B. Proof of Lemma \3\ 

Lower Bound: By expanding ( [Tol l around A = based on 
Taylor's series, it can be shown that 



p a L ut W = ^ L + o(x L+1 ) 



(20) 



where k is given in Lemma [3] 

Upper Bound: The asymptotic expansion of the upper 
bound in (fTTb is obtained as follows. To simplify notation, 
define B = dr a B~ x W . The term in ( fl2b can be written as 
P a (A) = E[A(P)] where the function A(6) is defined as 



A(6) = E{ft|G <6}Pr(G <b). 



where f2 = min 



C3 AG i 

(b-G-CiXG 2 - 6 ) 2 ' 



Mb) =E^n|G<Hpr(G<5 



The cumulative density function (CDF) of G is the probability 
that the number of nodes in the subset M. (g) is no more than 

L—l. Thus, Pr(G < g) = Zto ^iP" exp {-ciXg- s ) . 



Pr I - < G < b 



EiSl\G<l\px(G<± 



It can be expanded as 



E<! n\-<G <b} x 



Pr I - < G < b 



< 4c 3 A&~ 2 [1 + 0(A)] E 



G 2 - s \ G<- 



Pr [G < 



Ai(6) 



Pr(-<G<&], A^O. 



(21) 



A 2 (6) 

Next, by using the probability density function in (0]), 

(ciA)"- 1 r°° 



Mb) 



T(L) 



L-a 



Cl \(b/2)- s 



exp(— g)dg. (22) 



Note that for A — > 0, the integral in ( 1221 1 is bounded for L > a 
and unbounded for L < a. By separating these two cases, 



M 6 ) = (Cl w^ - r fi - a + 1, ciA (6/2)- 5 



Ai(6) < 



r(£) 

c L_l Ai _ 1(&/2) _ 5(L _ Q) 



r(L) 

Again, by using (0), 

! /-ciA(6/2) 



[1 + 0(A)], 



L > a 

L < a. 
(23) 



A 2 (6) 
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g L 1 cxp(-g)dg 
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M(&/2) 



-51(1,-1) 



r(L) 

2 5i (l - 2- 5 )c[ \ L b- SL 



g-CiAb-" _ e - Cl A(6/2) 



r(L) 

By combining (|24), (E2) and ED, 



E[A(5)] < < 



+ 0(A L+1 ) 



4E[B- 2 ]c 3 c^ i r (L - a + 1, c x A (6/2) 

<3(A Q+1 ), L>a 

_ 2 SL c[~Ei[B~ SL ] 



(24) 



-<5 



■A" + 



r(L) 



L < a 
where by using (O 



A L + 0(A L+1 ) 
(25) 
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The desired upper bound is obtained by combining ( fTTT i. ( |20b 
and 
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